The Kerr metric is an exact solution of the Einstein equations, generalizing the Schwarzschild metric to represent a spinning black hole. Chandrasekhar famously contradicted Landau’s assessment of the determination of this metric to claim that a derivation is “really very simple.” Surely only Chandrasekhar would consider his own presentation in the referenced volume “very simple”: I for one do not.

Knowing something of the features of the final metric, it is instructive to present a derivation that is simpler than Chandrasekhar’s and focuses on the mechanics of teasing out a solution from the Einstein equations. While the initial metric is every bit as general as that used by Chandrasekhar, the process will be guided by foreknowledge of the answer, with only evaluations made that lead to concrete progress. This is useful as an illustration of how to find the right thread to pull to unravel a rather complicated nonlinear problem.

Begin with a stationary axisymmetic metric of the form

$d{s}^2=(Adt+\epsilon Bd\phi {)}^2-(Cd\phi +\epsilon Ddt{)}^2-E^{2}d{r}^2-F^{2}d{\theta }^2$

where coefficients other than ε are functions of only the two spatial variables r and θ. Derivatives with respect to these variables will be denoted with subscripts 1 and 2, respectively.

The constant parameter $\epsilon$ measures deviation from a diagonal metric. In the limit $\epsilon \to 0$ the spherical symmetry of the Schwarzschild metric must be recovered:

$C\to rsin\theta F\to r,\epsilon \to 0$

In the same limit, A and E are functions of only the radial variable.

Choose basis forms in the Cartan formalism as

$\begin{array}{l}{\mathbf{\theta }}^0=Adt+\epsilon Bd\phi \\ {\mathbf{\theta }}^1=Edr\\ {\mathbf{\theta }}^2=Fd\theta \\ {\mathbf{\theta }}^3=Cd\phi +\epsilon Ddt\end{array}$

which are orthonormal by construction. The first and last basis forms are inverted with

$\begin{array}{c}dt=\frac{C}{G}{\mathbf{\theta }}^0-\epsilon \frac{B}{G}{\mathbf{\theta }}^3\\ d\phi =\frac{A}{G}{\mathbf{\theta }}^3-\epsilon \frac{D}{G}{\mathbf{\theta }}^0\end{array}G=AC-{\epsilon }^{2}BD$

The function combination G is the square root of the absolute value of the determinant of the submetric in the two variables t and φ.

Indices on forms are raised and lowered with the diagonal Minkowski metric [1,−1,−1,−1] .

The exterior derivatives of the basis forms are

$\begin{array}{l}d{\mathbf{\theta }}^0=A_{1}dr\wedge dt+A_{2}d\theta \wedge dt+\epsilon B_{1}dr\wedge d\phi +\epsilon B_{2}d\theta \wedge d\phi \\ \phantom{d{\mathbf{\theta }}^0}=-\frac{A_{1}C-{\epsilon }^2{B}_{1}D}{EG}{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^1-\frac{A_{2}C-{\epsilon }^2{B}_{2}D}{FG}{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^2\\ -\epsilon \frac{A_{1}B-B_{1}A}{EG}{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^3-\epsilon \frac{A_{2}B-B_{2}A}{FG}{\mathbf{\theta }}^2\wedge {\mathbf{\theta }}^3\\ \phantom{d{\mathbf{\theta }}^0}\equiv -P{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^1-R{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^2-H{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^3-K{\mathbf{\theta }}^2\wedge {\mathbf{\theta }}^3\\ \\ d{\mathbf{\theta }}^1=E_{2}d\theta \wedge dr=-\frac{E_2}{EF}{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^2\\ \\ d{\mathbf{\theta }}^2=F_{1}dr\wedge d\theta =\frac{F_1}{EF}{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^2\\ \\ d{\mathbf{\theta }}^3=C_{1}dr\wedge d\phi +C_{2}d\theta \wedge d\phi +\epsilon D_{1}dr\wedge dt+\epsilon D_{2}d\theta \wedge dt\\ \phantom{d{\mathbf{\theta }}^3}=\epsilon \frac{C_{1}D-D_{1}C}{EG}{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^1+\epsilon \frac{C_{2}D-D_{2}C}{FG}{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^2\\ +\frac{C_{1}A-{\epsilon }^2{D}_{1}B}{EG}{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^3+\frac{C_{2}A-{\epsilon }^2{D}_{2}B}{FG}{\mathbf{\theta }}^2\wedge {\mathbf{\theta }}^3\\ \phantom{d{\mathbf{\theta }}^3}\equiv J{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^1+L{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^2+Q{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^3+S{\mathbf{\theta }}^2\wedge {\mathbf{\theta }}^3\end{array}$

For convenience of evaluation, a comparative tabulation of the function combinations is useful:

$\begin{array}{l}H=\epsilon \frac{A_{1}B-B_{1}A}{EG}\\ J=\epsilon \frac{C_{1}D-D_{1}C}{EG}\\ K=\epsilon \frac{A_{2}B-B_{2}A}{FG}\\ L=\epsilon \frac{C_{2}D-D_{2}C}{FG}\end{array}\begin{array}{l}P=\frac{A_{1}C-{\epsilon }^2{B}_{1}D}{EG}\\ Q=\frac{C_{1}A-{\epsilon }^2{D}_{1}B}{EG}\\ R=\frac{A_{2}C-{\epsilon }^2{B}_{2}D}{FG}\\ S=\frac{C_{2}A-{\epsilon }^2{D}_{2}B}{FG}\end{array}$

These additional combinations look promising and will turn out to be important:

$P+Q=\frac{{\partial }_{1}lnG}{E}R+S=\frac{{\partial }_{2}lnG}{F}$

Consistent connection forms are identified with Cartan’s first equation of structure:

$d{\mathbf{\theta }}^i=-\mathbf{\omega }^i{}_k\wedge {\mathbf{\theta }}^k$

Attempting to do this by inspection tends to miss the crucial nonzero value for $\mathbf{\omega }^0{}_3$ . The straightforward procedure is to assume nonzero coefficients for all basis forms, expand the equation of structure and solve all nontrivial equations. The resulting connection forms are

$\begin{array}{l}\mathbf{\omega }^0{}_1=\mathbf{\omega }^1{}_0=P{\mathbf{\theta }}^0-\frac{H-J}{2}{\mathbf{\theta }}^3\\ \mathbf{\omega }^0{}_2=\mathbf{\omega }^2{}_0=R{\mathbf{\theta }}^0-\frac{K-L}{2}{\mathbf{\theta }}^3\\ \mathbf{\omega }^0{}_3=\mathbf{\omega }^3{}_0=\frac{H+J}{2}{\mathbf{\theta }}^1+\frac{K+L}{2}{\mathbf{\theta }}^2\\ \mathbf{\omega }^1{}_2=-\mathbf{\omega }^2{}_1=\frac{E_2}{EF}{\mathbf{\theta }}^1-\frac{F_1}{EF}{\mathbf{\theta }}^2\\ \mathbf{\omega }^1{}_3=-\mathbf{\omega }^3{}_1=-Q{\mathbf{\theta }}^3-\frac{H-J}{2}{\mathbf{\theta }}^0\\ \mathbf{\omega }^2{}_3=-\mathbf{\omega }^3{}_2=-S{\mathbf{\theta }}^3-\frac{K-L}{2}{\mathbf{\theta }}^0\end{array}$

It is simple to substitute these expressions back into the first equation of structure to verify the consistency of the set of forms.

It is now straightforward to evaluate curvature 2-forms using Cartan’s second equation of structure:

$\mathbf{R}^i{}_k=d\mathbf{\omega }_k^i+\mathbf{\omega }^i{}_m\wedge \mathbf{\omega }^m{}_k$

The curvature 2-forms are related to the Riemann tensor in so-called flat indices by

$\mathbf{R}^i{}_k=R^i{}_{klm}{\mathbf{\theta }}^l\wedge {\mathbf{\theta }}^m$

The Ricci tensor in the same indices is the contraction

$R_{ik}=R^m{}_{imk}$

Since vacuum solutions to the Einstein equation can be determined with either the Ricci or Einstein tensor, the former will be the goal of this presentation, bypassing listing all components of the Riemann tensor as unnecessary.

The six explicit equations for the curvature 2-forms are

$\begin{array}{l}\mathbf{R}^0{}_1=\mathbf{R}^1{}_0=d\mathbf{\omega }^0{}_1+\mathbf{\omega }^0{}_2\wedge \mathbf{\omega }^2{}_1+\mathbf{\omega }^0{}_3\wedge \mathbf{\omega }^3{}_1\\ \mathbf{R}^0{}_2=\mathbf{R}^2{}_0=d\mathbf{\omega }^0{}_2+\mathbf{\omega }^0{}_1\wedge \mathbf{\omega }^1{}_2+\mathbf{\omega }^0{}_3\wedge \mathbf{\omega }^3{}_2\\ \mathbf{R}^0{}_3=\mathbf{R}^3{}_0=d\mathbf{\omega }^0{}_3+\mathbf{\omega }^0{}_1\wedge \mathbf{\omega }^1{}_3+\mathbf{\omega }^0{}_2\wedge \mathbf{\omega }^2{}_3\\ \mathbf{R}^1{}_2=-\mathbf{R}^2{}_1=d\mathbf{\omega }^1{}_2+\mathbf{\omega }^1{}_0\wedge \mathbf{\omega }^0{}_2+\mathbf{\omega }^1{}_3\wedge \mathbf{\omega }^3{}_2\\ \mathbf{R}^1{}_3=-\mathbf{R}^3{}_1=d\mathbf{\omega }^1{}_3+\mathbf{\omega }^1{}_0\wedge \mathbf{\omega }^0{}_3+\mathbf{\omega }^1{}_2\wedge \mathbf{\omega }^2{}_3\\ \mathbf{R}^2{}_3=-\mathbf{R}^3{}_2=d\mathbf{\omega }^2{}_3+\mathbf{\omega }^2{}_0\wedge \mathbf{\omega }^0{}_3+\mathbf{\omega }^2{}_1\wedge \mathbf{\omega }^1{}_3\end{array}$

Running through the combinations of basis forms that can appear on the right-hand sides of these equations indicates that the third and fourth are simpler than the others. The possible resulting indices for the Riemann tensor can be represented as

$\begin{array}{l}R^0{}_{1[01,02,13,23]}=R^1{}_{0[01,02,13,23]}\\ R^0{}_{2[01,02,13,23]}=R^2{}_{0[01,02,13,23]}\\ R^0{}_{3[03,12]}=R^3{}_{0[03,12]}\\ R^1{}_{2[03,12]}=-R^2{}_{1[03,12]}\\ R^1{}_{3[01,02,13,23]}=-R^3{}_{1[01,02,13,23]}\\ R^2{}_{3[01,02,13,23]}=-R^3{}_{2[01,02,13,23]}\end{array}$

where alternatives are listed in the square subscripted brackets. This is already enough information to determine that four of the components of the Ricci tensor are identically zero,

$R_{01}=R_{02}=R_{13}=R_{23}=0$

since the components of the Riemann tensor that would compose them do not appear.

Now one must dig in and evaluate the curvature 2-forms one by one. Immediate simplifications using function combinations above are included: further simplifications are possible.

$\begin{array}{l}\mathbf{R}^0{}_1=\sum _{i=1,2}{\partial }_i(AP-\epsilon D\frac{H-J}{2})d{x}^i\wedge (\frac{C}{G}{\mathbf{\theta }}^0-\epsilon \frac{B}{G}{\mathbf{\theta }}^3)\\ +\sum _{i=1,2}{\partial }_i(\epsilon BP-C\frac{H-J}{2})d{x}^i\wedge (\frac{A}{G}{\mathbf{\theta }}^3-\epsilon \frac{D}{G}{\mathbf{\theta }}^0)\\ -(R{\mathbf{\theta }}^0-\frac{K-L}{2}{\mathbf{\theta }}^3)\wedge (\frac{E_2}{EF}{\mathbf{\theta }}^1-\frac{F_1}{EF}{\mathbf{\theta }}^2)\\ +(\frac{H+J}{2}{\mathbf{\theta }}^1+\frac{K+L}{2}{\mathbf{\theta }}^2)\wedge (Q{\mathbf{\theta }}^3+\frac{H-J}{2}{\mathbf{\theta }}^0)\\ \phantom{\mathbf{R}^0{}_1}=-[\frac{P_1}{E}+P^2+\frac{(H+3J)(H-J)}{4}+\frac{E_{2}R}{EF}]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^1\\ -[\frac{P_2}{F}+PR+\frac{(K+3L)(H-J)}{4}-\frac{F_{1}R}{EF}]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^2\\ -[\frac{H_1-J_1}{2E}+PH-QJ+\frac{E_2(K-L)}{2EF}]{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^3\\ -[\frac{H_2-J_2}{2F}+PK+\frac{S(H-J)-Q(K+L)}{2}-\frac{F_1(K-L)}{2EF}]{\mathbf{\theta }}^2\wedge {\mathbf{\theta }}^3\end{array}$

$\begin{array}{l}\mathbf{R}^0{}_2=\sum _{i=1,2}{\partial }_i(AR-\epsilon D\frac{K-L}{2})d{x}^i\wedge (\frac{C}{G}{\mathbf{\theta }}^0-\epsilon \frac{B}{G}{\mathbf{\theta }}^3)\\ +\sum _{i=1,2}{\partial }_i(\epsilon BR-C\frac{K-L}{2})d{x}^i\wedge (\frac{A}{G}{\mathbf{\theta }}^3-\epsilon \frac{D}{G}{\mathbf{\theta }}^0)\\ +(P{\mathbf{\theta }}^0-\frac{H-J}{2}{\mathbf{\theta }}^3)\wedge (\frac{E_2}{EF}{\mathbf{\theta }}^1-\frac{F_1}{EF}{\mathbf{\theta }}^2)\\ +(\frac{H+J}{2}{\mathbf{\theta }}^1+\frac{K+L}{2}{\mathbf{\theta }}^2)\wedge (S{\mathbf{\theta }}^3+\frac{K-L}{2}{\mathbf{\theta }}^0)\\ \phantom{\mathbf{R}^0{}_2}=-[\frac{R_1}{E}+PR+\frac{(H+3J)(K-L)}{4}-\frac{E_{2}P}{EF}]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^1\\ -[\frac{R_2}{F}+R^2+\frac{(K+3L)(K-L)}{4}+\frac{F_{1}P}{EF}]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^2\\ -[\frac{K_1-L_1}{2E}+RH+\frac{P(K-L)-S(H+J)}{2}-\frac{E_2(H-J)}{2EF}]{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^3\\ -[\frac{K_2-L_2}{2F}+RK-SL+\frac{F_1(H-J)}{2EF}]{\mathbf{\theta }}^2\wedge {\mathbf{\theta }}^3\end{array}$

$\begin{array}{l}\mathbf{R}^0{}_3={\partial }_2\left(E\frac{H+J}{2}\right)d\theta \wedge dr+{\partial }_1\left(F\frac{K+L}{2}\right)dr\wedge d\theta \\ -(P{\mathbf{\theta }}^0-\frac{H-J}{2}{\mathbf{\theta }}^3)\wedge (Q{\mathbf{\theta }}^3+\frac{H-J}{2}{\mathbf{\theta }}^0)\\ -(R{\mathbf{\theta }}^0-\frac{K-L}{2}{\mathbf{\theta }}^3)\wedge (S{\mathbf{\theta }}^3+\frac{K-L}{2}{\mathbf{\theta }}^0)\\ \phantom{\mathbf{R}^0{}_3}=-[PQ+RS+\frac{1}{4}(H-J{)}^2+\frac{1}{4}(K-L{)}^2]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^3\\ +\frac{1}{EF}\left[{\partial }_1\right(F\frac{K+L}{2})-{\partial }_2(E\frac{H+J}{2}\left)\right]{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^2\end{array}$

$\begin{array}{l}\mathbf{R}^1{}_2={\partial }_2\left(\frac{E_2}{F}\right)d\theta \wedge dr-{\partial }_1\left(\frac{F_1}{E}\right)dr\wedge d\theta \\ +(P{\mathbf{\theta }}^0-\frac{H-J}{2}{\mathbf{\theta }}^3)\wedge (R{\mathbf{\theta }}^0-\frac{K-L}{2}{\mathbf{\theta }}^3)\\ -(Q{\mathbf{\theta }}^3+\frac{H-J}{2}{\mathbf{\theta }}^0)\wedge (S{\mathbf{\theta }}^3+\frac{K-L}{2}{\mathbf{\theta }}^0)\\ \phantom{\mathbf{R}^1{}_2}=\frac{1}{2}\left[\right(Q-P\left)\right(K-L)+(R-S\left)\right(H-J\left)\right]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^3\\ -\frac{1}{EF}\left[{\partial }_1\right(\frac{F_1}{E})+{\partial }_2(\frac{E_2}{F}\left)\right]{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^2\end{array}$

$\begin{array}{l}\mathbf{R}^1{}_3=-\sum _{i=1,2}{\partial }_i(\epsilon DQ+A\frac{H-J}{2})d{x}^i\wedge (\frac{C}{G}{\mathbf{\theta }}^0-\epsilon \frac{B}{G}{\mathbf{\theta }}^3)\\ -\sum _{i=1,2}{\partial }_i(CQ+\epsilon B\frac{H-J}{2})d{x}^i\wedge (\frac{A}{G}{\mathbf{\theta }}^3-\epsilon \frac{D}{G}{\mathbf{\theta }}^0)\\ +(P{\mathbf{\theta }}^0-\frac{H-J}{2}{\mathbf{\theta }}^3)\wedge (\frac{H+J}{2}{\mathbf{\theta }}^1+\frac{K+L}{2}{\mathbf{\theta }}^2)\\ -(\frac{E_2}{EF}{\mathbf{\theta }}^1-\frac{F_1}{EF}{\mathbf{\theta }}^2)\wedge (S{\mathbf{\theta }}^3+\frac{K-L}{2}{\mathbf{\theta }}^0)\\ \phantom{\mathbf{R}^1{}_3}=[\frac{H_1-J_1}{2E}+PH-QJ+\frac{E_2(K-L)}{2EF}]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^1\\ +[\frac{H_2-J_2}{2F}-QL+\frac{R(H-J)+P(K+L)}{2}-\frac{F_1(K-L)}{2EF}]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^2\\ -[\frac{Q_1}{E}+Q^2-\frac{(3H+J)(H-J)}{4}+\frac{E_{2}S}{EF}]{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^3\\ -[\frac{Q_2}{F}+QS-\frac{(3K+L)(H-J)}{4}-\frac{F_{1}S}{EF}]{\mathbf{\theta }}^2\wedge {\mathbf{\theta }}^3\end{array}$

$\begin{array}{l}\mathbf{R}^2{}_3=-\sum _{i=1,2}{\partial }_i(\epsilon DS+A\frac{K-L}{2})d{x}^i\wedge (\frac{C}{G}{\mathbf{\theta }}^0-\epsilon \frac{B}{G}{\mathbf{\theta }}^3)\\ -\sum _{i=1,2}{\partial }_i(CS+\epsilon B\frac{K-L}{2})d{x}^i\wedge (\frac{A}{G}{\mathbf{\theta }}^3-\epsilon \frac{D}{G}{\mathbf{\theta }}^0)\\ +(R{\mathbf{\theta }}^0-\frac{K-L}{2}{\mathbf{\theta }}^3)\wedge (\frac{H+J}{2}{\mathbf{\theta }}^1+\frac{K+L}{2}{\mathbf{\theta }}^2)\\ +(\frac{E_2}{EF}{\mathbf{\theta }}^1-\frac{F_1}{EF}{\mathbf{\theta }}^2)\wedge (Q{\mathbf{\theta }}^3+\frac{H-J}{2}{\mathbf{\theta }}^0)\\ \phantom{\mathbf{R}^2{}_3}=[\frac{K_1-L_1}{2E}-SJ+\frac{P(K-L)+R(H+J)}{2}-\frac{E_2(H-J)}{2EF}]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^1\\ +[\frac{K_2-L_2}{2F}+RK-SL+\frac{F_1(H-J)}{2EF}]{\mathbf{\theta }}^0\wedge {\mathbf{\theta }}^2\\ -[\frac{S_1}{E}+QS-\frac{(3H+J)(K-L)}{4}-\frac{E_{2}Q}{EF}]{\mathbf{\theta }}^1\wedge {\mathbf{\theta }}^3\\ -[\frac{S_2}{F}+S^2-\frac{(3K+L)(K-L)}{4}+\frac{F_{1}Q}{EF}]{\mathbf{\theta }}^2\wedge {\mathbf{\theta }}^3\end{array}$

Some of the symmetries of the Riemann tensor pop out in these forms as written, but not all do. This can be used to advantage in determining relationships among metric functions.

Constructing the remaining components of the Ricci tensor is a simple matter of reading off components of the Riemann tensor. Components of the latter as written are antisymmetric in the first two indices when one is nonzero and always antisymmetric in the last two indices.

$\begin{array}{l}{R}_{00}=R^1{}_{010}+R^2{}_{020}+R^3{}_{030}\\ \phantom{R_{00}}=\frac{{\partial }_1\left(FP\right)}{EF}+\frac{{\partial }_2\left(ER\right)}{EF}+P(P+Q)+R(R+S)\\ +\frac{(H+J)(H-J)}{2}+\frac{(K+L)(K-L)}{2}\end{array}$

$\begin{array}{l}{R}_{11}=R^0{}_{101}+R^2{}_{121}+R^3{}_{131}\\ \phantom{R_{11}}=-\frac{P_1+Q_1}{E}-P^2-Q^2-\frac{E_2(R+S)}{EF}\\ -\frac{1}{EF}\left[{\partial }_1\right(\frac{F_1}{E})+{\partial }_2(\frac{E_2}{F}\left)\right]+\frac{(H-J{)}^2}{2}\end{array}$

$\begin{array}{l}{R}_{22}=R^0{}_{202}+R^1{}_{212}+R^3{}_{232}\\ \phantom{R_{22}}=-\frac{R_2+S_2}{F}-R^2-S^2-\frac{F_1(P+Q)}{EF}\\ -\frac{1}{EF}\left[{\partial }_1\right(\frac{F_1}{E})+{\partial }_2(\frac{E_2}{F}\left)\right]+\frac{(K-L{)}^2}{2}\end{array}$

$\begin{array}{l}{R}_{33}=R^0{}_{303}+R^1{}_{313}+R^2{}_{323}\\ \phantom{R_{33}}=-\frac{{\partial }_1\left(FQ\right)}{EF}-\frac{{\partial }_2\left(ES\right)}{EF}-Q(P+Q)-S(R+S)\\ +\frac{(H+J)(H-J)}{2}+\frac{(K+L)(K-L)}{2}\end{array}$

$\begin{array}{l}{R}_{03}=R^1{}_{013}+R^2{}_{023}=R^1{}_{310}+R^2{}_{320}=R_{30}\\ \phantom{R_{03}}=-\frac{{\partial }_1\left[F\right(H-J\left)\right]}{2EF}-\frac{{\partial }_2\left[E\right(K-L\left)\right]}{2EF}-PH+QJ-RK+SL\end{array}$

$\begin{array}{l}{R}_{12}=R^0{}_{102}+R^3{}_{132}=R^0{}_{201}+R^3{}_{231}=R_{21}\\ \phantom{R_{12}}=-\frac{P_2+Q_2}{F}-PR-QS+\frac{F_1(R+S)}{EF}+\frac{(H-J)(K-L)}{2}\\ \phantom{R_{12}}=-\frac{R_1+S_1}{E}-PR-QS+\frac{E_2(P+Q)}{EF}+\frac{(H-J)(K-L)}{2}\end{array}$

These last two expressions for the off-diagonal component are necessarily identical. Equating the apparently different parts leads to

$\begin{array}{c}{\partial }_2\left[E\right(P+Q\left)\right]={\partial }_1\left[F\right(R+S\left)\right]\\ {\partial }_2[{\partial }_{1}lnG]={\partial }_1[{\partial }_{2}lnG]\end{array}$

which in general is simply a statement of the commutativity of partial differentiation. While this necessarily holds and shows the consistency of the two expressions, it also suggests a path toward solution: if each side of this statement is separately zero, then the combination G of metric functions is separable in the two independent variables:

$G=AC-{\epsilon }^{2}BD=g\left(r\right)h\left(\theta \right)$

With this in mind, consider the combination of Ricci tensor components

$\begin{array}{l}{R}_{00}-R_{33}=\frac{{\partial }_1\left[F\right(P+Q\left)\right]}{EF}+\frac{{\partial }_2\left[E\right(R+S\left)\right]}{EF}\\ +(P+Q{)}^2+(R+S{)}^2\\ \phantom{R_{00}-R_{33}}=\frac{1}{EF}{\partial }_1(\frac{F}{E}{\partial }_{1}lnG)+\frac{1}{EF}{\partial }_2(\frac{E}{F}{\partial }_{2}lnG)\\ +\frac{1}{E^2}({\partial }_{1}lnG{)}^2+\frac{1}{F^2}({\partial }_{2}lnG{)}^2\\ R_{00}-R_{33}=\frac{1}{EF}{\partial }_1(\frac{F}{E}{\partial }_{1}lng)+\frac{1}{EF}{\partial }_2(\frac{E}{F}{\partial }_{2}lnh)\\ +\frac{1}{E^2}({\partial }_{1}lng{)}^2+\frac{1}{F^2}({\partial }_{2}lnh{)}^2\end{array}$

When equated to zero for a vacuum solution,

$\begin{array}{l}\frac{F}{E}{\partial }_1(\frac{F}{E}{\partial }_{1}lng)+\frac{F}{E}{\partial }_2(\frac{E}{F}{\partial }_{2}lnh)+(\frac{F}{E}{)}^2({\partial }_{1}lng{)}^2+({\partial }_{2}lnh{)}^2=0\end{array}$

one can achieve a separable equation by assuming

$\frac{F}{E}=f\left(r\right)$

Then the angular part of the equation, the second and fourth terms, can be set equal to a constant

${\partial }_2^{2}lnh+({\partial }_{2}lnh{)}^2=\frac{h_{22}}{h}=\mathrm{constant}$

and the function h is clearly a circular function in the compact variable θ. With the choice $h=sin\theta$ , the remainder of the equation is

$f{\partial }_1(f{\partial }_{1}lng)+f^2({\partial }_{1}lng{)}^2=\frac{f}{g}{\partial }_1\left(f{g}_1\right)=1$

The simplest way for this equation to hold is for the two functions to be equal, for which the solution is

$\begin{array}{c}{\partial }_1\left(f{f}_1\right)=\frac{1}{2}{\partial }_1^2{f}^2=1\\ f^2=g^2=r^2+\alpha r+\beta \equiv \Delta \end{array}$

and one can already begin to see the Kerr metric taking shape, just from assuming the simplest way toward a nontrivial solution.

Now look again at the statement of separability

$G=f\left(r\right)sin\theta =AC-{\epsilon }^{2}BD=AD(\frac{C}{D}-{\epsilon }^2\frac{B}{A})$

The difference of metric function ratios in brackets is unlikely to be equal to a separable product, but clearly is significant. Since it has to balance the product of two metric functions, denote it as a squared function:

${\rho }^2(r,\theta )=\frac{C}{D}-{\epsilon }^2\frac{B}{A}=c(r,\theta )-{\epsilon }^{2}b(r,\theta )$

The factor of the sine function can be assigned to both C and D in order to match the spherical symmetry of the Schwarzschild solution. Knowing also that the coefficient A in the Schwarzschild limit is a function of the radial variable, the factor f can be assigned with confidence to this metric function. With these assumptions, one has

$A=\frac{f}{\rho }B=\frac{bf}{\rho }C=\frac{csin\theta }{\rho }D=\frac{sin\theta }{\rho }$

where denominators have been assigned equally as the simplest choice possible.

With a focus on specific ratios of metric functions, the first set of intermediate metric function combinations becomes

$\begin{array}{l}H=-\epsilon \frac{A{\partial }_1(B/A)}{DE{\rho }^2}=-\epsilon \frac{b_1{f}^2}{{\rho }^{2}Fsin\theta }\\ J=\epsilon \frac{D{\partial }_1(C/D)}{AE{\rho }^2}=\epsilon \frac{c_{1}sin\theta }{{\rho }^{2}F}\\ K=-\epsilon \frac{A{\partial }_2(B/A)}{DF{\rho }^2}=-\epsilon \frac{b_{2}f}{{\rho }^{2}Fsin\theta }\\ L=\epsilon \frac{D{\partial }_2(C/D)}{AF{\rho }^2}=\epsilon \frac{c_{2}sin\theta }{f{\rho }^{2}F}\end{array}$

Differentiating the definition of ${\rho }^2$ ,

$2\rho {\rho }_1=c_1-{\epsilon }^2{b}_{1}2\rho {\rho }_2=c_2-{\epsilon }^2{b}_2$

the second set of metric function combinations can be written

$\begin{array}{l}P=\frac{f_1}{F}-\frac{(c_1+{\epsilon }^2{b}_1)f}{2{\rho }^{2}F}\\ Q=\frac{(c_1+{\epsilon }^2{b}_1)f}{2{\rho }^{2}F}\\ R=-\frac{(c_2+{\epsilon }^2{b}_2)}{2{\rho }^{2}F}\\ S=\frac{cot\theta }{F}+\frac{(c_2+{\epsilon }^2{b}_2)}{2{\rho }^{2}F}\end{array}$

All of the metric function combinations have a factor of F in the denominator. This means that when a vacuum equation is formed by for example setting $R_{03}$ equal to zero, then that particular metric function can be eliminated from the equation. This allows writing an equation that determines the interrelationship of b and c as functions of the two independent variables.

Such an equation cannot, unfortunately, be written in a completely integrable form. The difficulty essentially arises because the quantities

$\frac{c_1+{\epsilon }^2{b}_1}{{\rho }^2}=\frac{c_1+{\epsilon }^2{b}_1}{c-{\epsilon }^{2}b}\frac{c_2+{\epsilon }^2{b}_2}{{\rho }^2}=\frac{c_2+{\epsilon }^2{b}_2}{c-{\epsilon }^{2}b}$

are not in general integrable due to the sign difference between the numerator and denominator. One way to make these expressions integrable is to assume that the two function b and c each depend on only one of the independent variables, and there are two ways to do so:

${\rho }^2=c\left(\theta \right)-{\epsilon }^{2}b\left(r\right)\mathrm{or}{\rho }^2=c\left(r\right)-{\epsilon }^{2}b\left(\theta \right)$

For the first choice, $J=K=0$.  The vacuum equation formed by setting $R_{03}$ equal to zero is in that case

$\begin{array}{r}-\frac{{\partial }_1\left(FH\right)}{2EF}+\frac{{\partial }_2\left(EL\right)}{2EF}-PH+SL=0\\ \frac{f}{2}{\partial }_1\left(\frac{b_1{f}^2}{{\rho }^{2}sin\theta }\right)+\frac{f}{2}{\partial }_2\left(\frac{c_{2}sin\theta }{f^2{\rho }^2}\right)\hfill \\ +\frac{b_1{f}^2}{{\rho }^{2}sin\theta }(f_1-\frac{{\epsilon }^2{b}_{1}f}{2{\rho }^2})+\frac{c_{2}sin\theta }{f{\rho }^2}(cot\theta +\frac{c_2}{2{\rho }^2})=0\\ {\partial }_1\left(b_1{f}^4\right)+\frac{1}{sin\theta }{\partial }_2\left(c_2{sin}^3\theta \right)=0\end{array}$

which is a separable equation. Setting the second term equal to a constant implies

${\partial }_2\left(c_2{sin}^3\theta \right)\sim sin\theta \to c\left(\theta \right)\sim \frac{1}{{sin}^2\theta }$

This solution would not match the Schwarzschild solution in the limit $\epsilon \to 0$ and so cannot be correct. Turn to the second choice, for which $H=L=0$.  The vacuum equation formed by setting $R_{03}$ equal to zero is in this case

$\begin{array}{r}\frac{{\partial }_1\left(FJ\right)}{2EF}-\frac{{\partial }_2\left(EK\right)}{2EF}+QJ-RK=0\\ \frac{f}{2}{\partial }_1\left(\frac{c_{1}sin\theta }{{\rho }^2}\right)+\frac{f}{2}{\partial }_2\left(\frac{b_2}{{\rho }^{2}sin\theta }\right)\\ +\frac{c_{1}sin\theta }{{\rho }^2}\frac{c_{1}f}{2{\rho }^2}-\frac{b_{2}f}{{\rho }^{2}sin\theta }\frac{{\epsilon }^2{b}_2}{2{\rho }^2}=0\\ {\partial }_1^{2}c+\frac{1}{sin\theta }{\partial }_2\left(\frac{b_2}{sin\theta }\right)=0\end{array}$

This equation is also separable and simpler than that of the first choice. The function b can written either in terms of sine or cosine. Selecting the former one has

$b\left(\theta \right)={sin}^2\theta c\left(r\right)=r^2+\gamma r+\delta$

To determine the constants of integration, an independent vacuum equation with F eliminated is needed. Consider the combination of Ricci tensor components

$\begin{array}{l}{R}_{00}+R_{33}=\frac{{\partial }_1\left[F\right(P-Q\left)\right]}{EF}+\frac{{\partial }_2\left[E\right(R-S\left)\right]}{EF}\\ +(P-Q)(P+Q)+(R-S)(R+S)\\ +(H-J)(H+J)+(K-L)(K+L)\\ \phantom{R_{00}+R_{33}}=\frac{f{\partial }_1\left[F\right(P-Q\left)\right]}{F^2}+\frac{{\partial }_2\left[F\right(R-S\left)\right]}{F^2}\\ +(P-Q)\frac{f{\partial }_{1}lnf}{F}+(R-S)\frac{{\partial }_{2}lnsin\theta }{F}-J^2+K^2\\ R_{00}+R_{33}=\frac{{\partial }_1\left[fF\right(P-Q\left)\right]}{F^2}+\frac{{\partial }_2\left[F\right(R-S)sin\theta ]}{F^{2}sin\theta }-J^2+K^2\end{array}$

Setting this equal to zero and substituting the angular function b gives

$\begin{array}{r}{\partial }_1[f{f}_1-\frac{c_1{f}^2}{{\rho }^2}]-\frac{1}{sin\theta }{\partial }_2[cos\theta +\frac{{\epsilon }^2{b}_2}{{\rho }^2}sin\theta ]-\frac{{\epsilon }^2{c}_1^2{sin}^2\theta }{{\rho }^4}+\frac{{\epsilon }^2{b}_2^2{f}^2}{{\rho }^4{sin}^2\theta }=0\\ 2{\rho }^4-{\rho }^4{\partial }_1\left(\frac{c_1{f}^2}{{\rho }^2}\right)-\frac{2{\epsilon }^2{\rho }^4}{sin\theta }{\partial }_2\left(\frac{{sin}^2\theta cos\theta }{{\rho }^2}\right)-{\epsilon }^2{c}_1^2{sin}^2\theta +4{\epsilon }^2{f}^2{cos}^2\theta =0\end{array}$

This looks messy, but by grouping on powers of sine it simplifies somewhat to

$\begin{array}{l}\left[{\partial }_1\right(c_1{f}^2)-c_1^2-4f^2+2c]{\epsilon }^2{sin}^2\theta \\ +4{\epsilon }^2(f^2-c)-c{\partial }_1\left(c_1{f}^2\right)+c_1^2{f}^2+2c^2=0\end{array}$

While further simplification is possible here, it is straightforward to substitute the radial functions f and c and do a bit of algebra. There are significant cancellations in the process:

$\begin{array}{l}\left[\gamma \right(\alpha -\gamma )-2(\beta -\delta \left)\right]{\epsilon }^2{sin}^2\theta \\ +\left[2\right(\beta -\delta )-\gamma (\alpha -\gamma \left)\right]r^2\\ +\left[4\alpha \right({\epsilon }^2-\delta )+2\gamma (\beta +\delta -2{\epsilon }^2\left)\right]r\\ +\left[2\right(\beta -\delta \left)\right(2{\epsilon }^2-\delta )-\gamma (\alpha \delta -\beta \gamma \left)\right]=0\end{array}$

Each of the quantities in square brackets must be separately zero for this to hold. Two of the resulting equations are identical, so only three of the four constants can be determined. From the given grouping it is clear that a good choice is $\beta =\delta ={\epsilon }^2$ . For a third constant one can choose either $\gamma =0$ or $\gamma =\alpha$ : the Kerr metric uses the former, but clearly a slightly different solution is possible.

For the standard Kerr metric, the two radial functions are now

$f\left(r\right)=r^2+\alpha r+{\epsilon }^{2}c\left(r\right)=r^2+{\epsilon }^2$

All that remains is to determine the form of the metric function F relative to the others. For that one could consider the combination of Ricci tensor components

$\begin{array}{l}-R_{11}-R_{22}-R_{00}+R_{33}\\ =\frac{1}{EF}\left[{\partial }_1\right(\frac{F_1}{E})+{\partial }_2(\frac{E_2}{F}\left)\right]-2PQ-2RS-\frac{(H-J{)}^2}{2}-\frac{(K-L{)}^2}{2}\end{array}$

but the expression in square brackets is less simple than it appears. Instead use a sum of the two forms of the off-diagonal component:

$\begin{array}{l}-2R_{12}=\frac{F}{E}{\partial }_1\left(\frac{R+S}{F}\right)+\frac{E}{F}{\partial }_2\left(\frac{P+Q}{E}\right)+2PR+2QS-(H-J)(K-L)\\ \phantom{-2R_{12}}=fcot\theta {\partial }_1\left(\frac{1}{F^2}\right)+f_1{\partial }_2\left(\frac{1}{F^2}\right)-\frac{{\epsilon }^2{b}_2{f}_1}{{\rho }^2{F}^2}+\frac{c_{1}fcot\theta }{{\rho }^2{F}^2}\end{array}$

The somewhat surprising reduction in the number of terms here is because the sine factors in the nonzero functions J and K cancel against each other. This simplification would not occur if H and L were also nonzero.

Equating to zero for a vacuum solution and regrouping,

$fcot\theta \left[{\partial }_1\right(\frac{1}{F^2})+\frac{c_1}{{\rho }^2{F}^2}]+f_1\left[{\partial }_2\right(\frac{1}{F^2})-\frac{{\epsilon }^2{b}_2}{{\rho }^2{F}^2}]=0$

it is clear by inspection that F and ρ are identical functions, which completes the determination of all metric functions. A final check of consistency would be to evaluate all components of the Ricci tensor individually to confirm that they are all zero. This can be done with a computer algebra system for efficiency and will not be included here.

To summarize, the coefficient functions determined are

$\begin{array}{l}A=\frac{f}{\rho }B=\frac{f{sin}^2\theta }{\rho }C=\frac{(r^2+{\epsilon }^2)sin\theta }{\rho }D=\frac{sin\theta }{\rho }\\ \\ E=\frac{\rho }{f}F=\rho \\ \\ {\rho }^2=r^2+{\epsilon }^2{cos}^2\theta f^2=r^2+\alpha r+{\epsilon }^2=\Delta \end{array}$

The physical interpretation of the two constants can be achieved by comparing an expansion of the metric in inverse powers of r to the known weak-field limit of a rotating body, given for example as Equation (2) of this excerpt. The coefficient of $d{t}^2$ as developed here is

$A^2-{\epsilon }^2{D}^2=\frac{r^2+\alpha r+{\epsilon }^2{cos}^2\theta }{{\rho }^2}=1+\frac{\alpha r}{{\rho }^2}=1+\frac{\alpha }{r}+O\left(\frac{1}{r^2}\right)$

The numerator of the second term must be equal to the negative of twice the gravitational mass of the body: $\alpha =-2m$ . The coefficient of the cross term $dtd\phi$ as developed here is

$2\epsilon (AB-CD)=\frac{2\alpha \epsilon r{sin}^2\theta }{{\rho }^2}=\frac{2\alpha \epsilon {sin}^2\theta }{r}+O\left(\frac{1}{r^2}\right)$

The constant part of the numerator here must equal four times the angular momentum of the body:

$2\alpha \epsilon =4J\to \epsilon =-\frac{J}{m}=-a$

With these identifications, the metric functions become

$\begin{array}{l}A=\frac{f}{\rho }B=\frac{f{sin}^2\theta }{\rho }C=\frac{(r^2+a^2)sin\theta }{\rho }D=\frac{sin\theta }{\rho }\\ \\ E=\frac{\rho }{f}F=\rho \\ \\ {\rho }^2=r^2+a^2{cos}^2\theta f^2=r^2-2mr+a^2=\Delta \end{array}$

The complete Kerr metric with these coefficient functions is

$d{s}^2=\frac{\Delta }{{\rho }^2}[dt-a{sin}^2\theta d\phi {]}^2-\frac{{sin}^2\theta }{{\rho }^2}\left[\right(r^2+a^2)d\phi -adt{]}^2-\frac{{\rho }^2}{\Delta }d{r}^2-{\rho }^{2}d{\theta }^2$

Seeing the metric expressed this way in this discussion is what inspired the entire derivation. It struck me that such evenly matched powers of ${\rho }^2$ could only balance one another through ratios of metric functions. From there it was simple to assume similar behavior for the functions leading to cross terms, so that paying attention to ratios was the key to the entire process.

As stated at the outset, the initial metric form is every bit as general as that employed by Chandrasekhar. The inclusion of two cross terms instead of one makes the whole process as well as the final statement of the metric much simpler, in my opinion. Not that it is particularly simple...

Uploaded 2018.02.10 — Updated 2018.04.06 analyticphysics.com